Velocity Dependence of Quarkonium Dissociation Temperature in High-Energy 

Nuclear Collisions 
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In high-energy nuclear collisions, heavy quark potential at finite temperature controls the quarko- 
nium production. Including the relaxation of the medium induced by the relative velocity between 
quarkonia and the deconfined expanding matter, the Debye screening is reduced and the quarko- 
nium dissociation takes place at a higher temperature. As a consequence of the velocity dependent 
dissociation temperature, the quarkonium suppression at high transverse momentum is significantly 
weakened in high energy nuclear collisions at RHIC and LHC. 



Heavy quarkonia J/ip and T are tightly bound 
hadronic states. Their dissociation temperature Td is, in 
general, higher than the critical temperature T c for the 
deconfinement phase transition [l[ in high-energy nuclear 
collisions 0-01 . Therefore, the measured cross sections 
of quarkonia carry the information of the early stage hot 
and dense medium. They have long been considered as 
a signature of the formation of the new state of matter, 
the so-called quark-gluon plasma si. 

The quarkonium dissociation in a static deconfined 
quark matter is generally described in terms of the 
screening effect. The heavy quark potential, which is 
normally taken as the Cornell form [7[ and can be cal- 
culated through a non-relativistic quantum chromody- 
namic potential [|| and lattice simulations @ , is reduced 
to a Yukawa-like potential due to the Debye screening. 
When the screening radius becomes smaller than the 
quarkonium size, the bound state dissociates. Substitut- 
ing the screened potential, extracted from lattice simula- 
tions l3, U | , into the Schrodinger equation for the wave 
function of the quarkonium state, one obtains the disso- 
ciation temperature that corresponds to the zero hinging 
energy and infinite size of the di-quark system 12|, 1 1 31 ] . 
For charmonia, while the excited states \c and tp' start 
to dissociate already around T c , the calculated dissocia- 
tion temperature for the ground state J/tjj is much higher 
than the critical temperature 12 1. 

The quarkonia produced in relativistic heavy ion col- 
lisions are, however, not at rest in the medium. There 
exists a relative velocity between the quarkonia and the 
expanding medium. The question is what is the velocity 
dependence of the heavy quark potential at finite tem- 
perature. The screening effect is due to the rearrange- 
ment of the charged particles when a pair of heavy quarks 
(source) presents in the medium. For a moving source, 
it will take a longer time for the source to interact with 
the medium, comparing with that of a stationary source. 
This 'delay' of the response reduces the screening charges 
around the source and thus weakens the screening effect. 
In relativistic heavy ion collisions, the averaged trans- 



verse momentum of the initially produced J/ifts is about 
several GeV (T^Ji corresponding to an averaged relative 
velocity of the order of (0.5 — l)c. A significant modifi- 
cation of the Debye screening is expected for such fast 
moving J/ips, especially for tho se p roduced in higher 
transverse momentum region 1^, 16 1. In this Letter, we 
study the velocity dependence of the heavy quark po- 
tential and the quarkonium dissociation temperature in 
a transport approach. The velocity induced effects on 
charmonium production at both Relativistic Heavy Ion 
Collider (RHIC) and Large Hadron Collider (LHC) will 
be discussed. In the following calculations we take the 
speed of light c = 1. 

For a static source located at r = 0, the ambient charge 
density po( r ) is modified by the screening potential Vq(t) 
at finite temperature T |17| . 
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- qi V (T)/T 
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where qi is the charge of the particles of the specie i, 
and fi is the initial particle density without the source. 
The neutrality condition for the total charge has been 
considered here. 

For a moving source in the medium with relative ve- 
locity v, the non-equilibrium charge density p(r, t) in the 
source-rest frame satisfies the transport equation in the 
relaxation time approximation, 



d t p - v • Vp = - {p - pa) /t, 



(2) 



where r is the relaxation time of the medium. Tak- 
ing the limit t — > oo, the final distribution pf(r, L) = 
limt^oo p(r, t) becomes stable and is characterized by the 
equation 



L ■ Vp f = p f - pa, 



(3) 



where we have introduced the relaxation length defined as 
L = vr which controls the velocity dependence of the De- 
bye screening. Since the screening charge distribution is 
proportional to the screening potential, see equation (fJJ , 
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the potentials Vq and V corresponding to a stationary 
and moving source, respectively, satisfy the same equa- 
tion ([3]) which can be solved analytically, 



/>oo 

V(r, L) = / V Q (r + \L)t 
Jo 



" x d\. 



(4) 



It is obvious that for a static source with L = we have 
V(v,0) = Vo(r). 

With the potentials Vq and V, the screening radius r d 
can be expressed as 

1 J rf 3 r r P/ (r, L) 1 J ri 3 r r7(r, L) 
rdl j 2/d 3 r P/ (r,L) 2 /d3rF (r) ' W 

where the total charge conservation has been accounted 
for. For a spherically symmetric potential Voir), the inte- 
grations in ([5]) can be greatly simplified and the effective 
radial potential at nonzero velocity can be factorized as 



V(r,L) = V (r)W(r/L) 



with the modification factor 



W(y) 



1 



2 + y 3 Z(y) - (y 2 - y + 2)e~« 
3y 2 



(6) 



(7) 



where the function Z(y) is defined as Z(y) = dt e~ t /t. 

Now we apply the above transport solutions to the 
quarkonium dissociation in hot and dense matter cre- 
ated in high-energy nuclear collisions. The interaction 
between two quarks in vacuum can be well characterized 
by the Cornell potential [l2[ Vo(r) = —a/r + or with 
coupling constant a = it/ 12 and string tension a = 0.2 
GeV 2 . At finite temperature, the screening potential 
for a stationary pair of heavy quarks can be written 
as El El 



Voir) 



-fir 



2*r(f) 



1/2 



^i((H 2 ), (8) 



where T and K are the Gamma and modified Bessel func- 
tions. The temperature of the medium is hidden in the 
screening mass /i(T) which can be extracted fl2l | from 
lattice QCD calculated free energy 11 1. 



To establish a unique mapping between the relaxation 
length and the velocity, we estimate the relaxation time 
of the hot and dense matter by considering its electric 
analogue. When an electric charge is put into a conduct- 
ing medium, the charge is neutralized in a time scale of 
r = 1/ (4:ir<T e a e ), where a e is the electric conductivity of 
the medium, and a e is the fine-structure constant. We re- 
place a* by the color conductivity <r c rj 0.4T, estimated 
from [l9j], and a e by a, the relaxation length becomes 
L = l5v/(2ir 2 T). 

In FigQ] one sees the velocity induced change in the 
heavy quark potential at a fixed temperature T = 1.5T C . 
The stationary potential is taken from the lattice simula- 
tion © . Since the screening length is proportional to the 



velocity and inversely proportional to the temperature of 
the medium, the potential well becomes deeper and the 
screening effect becomes weaker, when the quarkonium 
velocity relative to the medium increases. 
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FIG. 1. The velocity dependence of the screening potential 
at a fixed temperature T = 1.5T C . a is the string tension and 
the stationary potential is taken as the free energy [Tol - H^ ]. 

With the known potentials Vo(r) and V(r,L), the 
screen radius Vd{T, L) at finite temperature can be calcu- 
lated through ([5]). In the rest frame of the di-quark sys- 
tem, the condition for dissociating a quarkonium should 
not depend on its relative velocity, namely the critical 
screening radius is a constant, 



r d {T d {L),L) 



1 1 

r d (T d (0),0) = 



161^(3/4) op 
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(9) 



with /i(T d (0)), which determines the relationship between 
the dissociation temperatures T d (L) for a moving quarko- 
nium and T d (0) for a static quarkonium. When T^O) 
runs from T c to 2.5T C , the screening radius r d (T d (0),0) 
runs from 0.44 fm to 0.21 fm. The running dissociation 
temperature T d (v) is shown in Fig|2l Once we fix T d (0), 
its velocity dependence can be obtained from the corre- 
sponding curve. As expected, when a quarkonium moves 
at a large velocity relative to the medium, the screening 
effect becomes weaker and the dissociation temperature 
becomes higher. The velocity induced shift of the disso- 
ciation temperature can be as large as AT d (v)/T c ~ 1 for 
fast quarkonia. For charmonium, the dissociation tem- 
perature is T d /T c ~ 1 — 2atw = 12, EH but goes up 
to 1.2 — 2.7 at v = 0.8, see Fig|2] Considering the fact 
that the fireball temperature formed in heavy ion colli- 
sions at RHIC energy is in the region T/T c ~ 1 — 2, the 
J/ip transverse momentum spectrum should be sensitive 
to the velocity dependence of the dissociation tempera- 
ture. The much higher dissociation temperature for fast 
moving J ftps will lead to a weaker suppression in the high 
transverse momentum region. 

In order to quantitatively see the effect of the running 
dissociation temperature T d (v) on quarkonium produc- 
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FIG. 2. The dissociation temperature Td(v), starting at dif- 
ferent stationary values 1 < T<j(0)/T c < 2.5, T c is the critical 
temperature of the quark matter. 



tion in high-energy nuclear collisions, we take a detailed 
transport approach [20l ] to describe the dynamical evo- 
lution of the hot and dense medium. The model con- 
tains transport equations for the quarkonium motion in 
the medium and hydrodynamic equations for the space- 
time evolution of the medium. In heavy ion collisions 
there are two sources for quarkonium production: the 
primordial production at the initial state and the regen- 
eration in the hot medium. During the evolution, all 
of the produced quarkonia suffer from the medium in- 
duced dissociation, dominantly by the gluon interactions. 
The model used here describes well both J/ip [2l|, 22 1 
and T 



23( productions. In order to demonstrate the 



velocity effect on the quarkonia production, it is neces- 
sary to study the transverse momentum distributions. 
We consider the differential nuclear modification fac- 
tor RAA(pt) = NAA(pt)/ (NcoiiNpp{pt)) as a function 
of quarkonium transverse momentum p t . NAA(pt) and 
N pp (pt) are differential quarkonium yields in heavy ion 
and elementary p+p collisions, and N co u is the number 
of nucleon+nucleon collisions in heavy ion collisions. We 
will focus on the high p t behavior of RAA(pt)- 

Fig El shows the J/ip Raa{pi) for a constant and a 
running dissociation temperature in central Au+Au col- 
lisions at top RHIC energy ^Jsnn = 200 GeV. As one 
can see, at the low p t region (< 3.5 GeV), the experi- 
mental results of J/ip RAA(pt) is less than 0.4. However, 
in the higher p t region, the value of the nuclear modifi- 
cation factor becomes higher RAA{pt) ~ 0.6 indicating 
weaker suppression in J/ip production. Note that this pt 
dependence can not be reproduced by a constant disso- 
ciation temperature unless a strong Cronin effect even at 
extremely high p t (26|. As discussed in j27|, the strong 
Cronin effect at high p t region is not favored by the latest 
J/ip data from d+Au collisions. In our calculation, the 
Cronin effect has been characterized by a Gauss smearing 
scheme [28| , which contribute little to the high p t region 



Pt > 6 GeV. For the stationary charmonia, the dissocia- 
tion temperature calculated from the Schrodingcr equa- 
tion is in between (l.f — 2.I)T C for the ground state J/ip 
and T c for the excited states ip' an d Y c , depending on 
the used heavy quark potential [l2. 29l 3fjj. 

Since the fireball temperature in a central collision is 
much higher than T c , almost all the excited states are 
dissociated in the medium, we will consider mainly the 
ground state. Considering the fact that the contribution 
from the decay of the excited states to final J/ ip is about 
40%, there is an up limit of 60% for J/ip R A a- In FigJ3] 
the three dashed lines represent the results with a con- 
stant dissociation temperature Td =1.3, 1.6, and 1.9T C 
from bottom to top, respectively. As one can see in the 
figure, the numerical results with a constant T c i are all 
much less than 0.6 and underestimate the J/ip Raa- 
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FIG. 3. The J/ip nuclear modification factor RAA{pt) at 
RHIC. The data are from the PHENIX @] and STAR [H 
collaborations, the solid line is the calculation with a running 
dissociation temperature starting at Td(0)/T c — 1.6, and the 
dashed lines are the calculations with a constant dissociation 
temperature. 

We now analyze the results with the running disso- 
ciation temperature, see the solid line in Fig(3] To fit 
the J/ip yield, we take T d (0) = 1.6T C . Note that for an 
expanding fireball, p t = in the laboratory frame corre- 
sponds generally to a nonzero velocity in the rest frame 
of the fireball, therefore the running dissociation temper- 
ature even at p t = is already affected by the velocity 
v. That is why the Raa at p t — does not coincide 
with the calculation with a constant Td = 1.6T C . From 
FigHJ one sees that the running dissociation tempera- 
ture quickly increases as a function of the velocity. On 
the other hand, the maximum temperature of the fire- 
ball in a central Au+Au collisions is T max ps 2T C 31 , 32| . 



Therefore, at sufficient high transverse momentum, the 
dissociation temperature may stay above the maximum 
temperature. As a result, those high p t J ftps will sur- 
vive in the quark gluon plasma. For the same reason, 
those high p t excited states ip's and XcS produced in the 
more peripheral region where the temperature is lower 
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than their dissociation temperature will also survive. The 
competition between the velocity (pt) dependent dissoci- 
ation temperature and the fireball temperature leads to 
Raa > 60% at high p t , as shown in Fig|3J It is clear 
in the figure that our calculation with the running disso- 
ciation temperature is consistent with the experimental 
observation. High statistics data are needed in order to 
confirm this ansatz for quarkonium productions in high- 
energy nuclear collisions. 

In order to further test the model, the p t dependence of 
the RAA{pt) for prompt J /ips from minimum bias Pb+Pb 
collisions at LHC energy ^/sTTn = 2.76 TeV is also cal- 
culated in our model. The results are compared with 
the experimental data 33[ in FigSJ For this calcula- 
tion, the charm quark production cross section is taken 
as daNN/dy 



Sig- 



0.62 mb at middle rapidity [34 
nificant regeneration has been reported with the large 
charm cross section 3J, [35j . Since the charm quarks in- 
teract strongly with the medium, losing its initial energy, 
the regenerated charmonia are soft, leading to a large 
Raa in the low p t region. Similar to the case at RHIC, 
the CMS experimental value of Raa ~ 0.3 at high p t 
can be reproduced only when the running dissociation 
temperature is considered, while all other results with a 
constant underpredict the values of Raa • At the same 
Pt ~ 6 GeV region, the value of Raa from LHC is much 
lower than that from RHIC, implying that a much hotter 
medium has been formed in heavy ion collisions at the 
higher energy. 




FIG. 4. The J/ip nuclear modification factor RAA(pt) at 
LHC. The data are from the CMS Collaboration [23], the 
solid line is the calculation with a running dissociation tem- 
perature starting at Ta(0)/T c — 1.6, and the dashed lines are 
the calculations with a constant dissociation temperature. 

In summary, we studied the heavy quark potential and 
dissociation temperature for moving quarkonia in quark 
gluon plasma in high-energy nuclear collisions. For a 
moving heavy quark pair in the hot medium, the screen- 
ing potential is reduced and the dissociation temperature 
is enhanced. As a consequence of the running dissoci- 
ation temperature in charmonium production, the J/tj) 



suppression becomes significantly weaker at high trans- 
verse momentum. 
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